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We propose a torsional response raised by lattice dislocation in Weyl semimetals akin to chiral 
magnetic effect; i.e. a fictitious magnetic field arising from screw or edge dislocation induces charge 
current. We demonstrate that, in sharp contrast to the usual chiral magnetic effect which vanishes 
in real solid state materials, the torsional chiral magnetic effect exists even for realistic lattice 
models, which implies the experimental detection of the effect via SQUID or nonlocal resistivity 
measurements in Weyl semimetal materials. 
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Recently, many candidate materials for Dirac semimet¬ 
als and Weyl semimetals (WSMs) [1-4], have been discov¬ 
ered [5-23]. These topological semimetals are intriguing 
because of exotic transport phenomena associated with 
the chiral anomaly in quantum field theory [24], such 
as the anomalous Hall effect [25, 26], chiral magnetic ef¬ 
fect (CME) [27], negative longitudinal magnetoresistance 
[5, 12, 19-21, 28, 29], and chiral gauge field [30]. 

Among them, the CME has been discussed in broad 
areas of quantum many-body physics, including nuclear 
and nonequilibrium physics as well as condensed matter 
physics. It is the generation of charge current parallel to 
an applied magnetic field even in the absence of electric 
fields. In nuclear physics, together with the chiral vor¬ 
tical effect [31], it is expected to play an important role 
in heavy ion collisions experiments [32, 33]. The CME 
also caused a stir in nonequilibrium statistical physics, 
since it leads to the existence of the ground state which, 
recently, attracts a renewed interest in connection with 
the realization of quantum time crystal [34], and then the 
CME has been studied from this point of view [35, 36]. 
However, unfortunately, their results are negative for its 
realization: the macroscopic ground state current in re¬ 
alistic WSMs is always absent. 

In this letter, we propose a chiral response in WSMs, 
named “torsional chiral magnetic effect (TCME)”, in 
which the ground state charge current is caused by the 
effective magnetic field induced by lattice dislocation as 
shown in EIG.l. By using the Cartan formalism of the 
differential geometry, we can describe the lattice strain 
and dislocation in terms of vielbein and torsion [37]. 
Erom the viewpoint of the quantum field theory in curved 
space-time, the TCME is raised by the mixed action of 
electromagnetic and torsional fields that is prohibited in 
four-dimensional spacetime with the Lorentz symmetry, 
but made possible in non-relativistic band electrons in 
solid state systems. Eurthermore, we demonstrate that 
the TCME is possible in realistic lattice models by carry¬ 
ing out numerical calculations. Our results imply the ex¬ 
istence of experimentally observable current induced by 
the TCME in real WSM materials. We also resolve the 


relation between our results and the no-go theorem that 
the CME is absent in equilibrium states [35, 36]. Eirst of 
all, we clarify the notations. The indices = x, y, z 

and a,b, - ■ ■ = x,y,z represent the coordinates in the lab¬ 
oratory and local orthogonal (or Lorentz) frames, respec¬ 
tively. In the following, we use the Einstein summation 
convention. 

Linear response theory for torsional response — Here, 
we briefly introduce the Cartan formalism, which can 
be applied to description of crystal systems with lattice 
strain as follows. It is an approach to curved space and 
based on the local orthonormal frame form e“ = ef{r)dr^, 
where the coefficient fields e“(r) are referred to as the 
vielbein [38]. We introduce the coordinate measured by 
an observer on the deformed lattice and the labora¬ 
tory coordinate r*, and identify its exterior derivative 
as the local orthonormal frame e“ = Then, to 

the first order in the displacement field u, the vielbein 
is written as e“ = bf — du°‘jdK. Eor this observer, 
the lattice is not deformed, and then the Hamiltonian 
of the system is given by H[—\dRx,—idiiv,—idRs) = 
H{—ie^dri,—icyd^i,—ie\dri), where H{px,py,pz) is the 
Hamiltonian without lattice deformation and e^ is the 
inverse of e“. In this way, the emergent vielbein appears, 
and therefore we can describe the elastic response by us¬ 
ing the Cartan formalism. The coupling between the 
vielbein and electrons is similar to the minimal coupling 
of the U{1) gauge held, Pi ^ Pi — eAi. Then, we can de- 
hne the analog of the held strength by T“- = dicf — djef, 
which is referred to as the torsion, or “torsional mag¬ 
netic held” (TMF) [39-41], where the spin connection is 
dropped for simplicity. Using the displacement vector, 
the torsion is rewritten as r“ = {djdi — didj)u°‘. The 
point is that, if u“(r) is a well-dehned function, the tor¬ 
sion is always zero, and the multivaluedness of u“(r) is 
necessary for nonzero torsion. Indeed, the edge dislo¬ 
cation along z-axis with Burgers vector bgjf causes the 
TMF, T^y = —bgS^‘^\x,y), and the screw one with bgZ, 
Tfy = -bgS^‘^'>{x,y), as shown in FIG.l. For more details 
about the lattice strain and differential geometry, see, for 
example, Refs. [42-44]. 
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Now, using the linear response theory with the Car- 
tan formalism, we investigate the TCME of WSMs due 
to dislocation. We calculate the current density in the 
presence of TMF and magnetic field up to the linear or¬ 
der. We use the model of a pair of Weyl fermions with 
the opposite chirality, whose Weyl points are at k = \^ 
and in the momentum space, and Fermi energies are 
given by E = vf^q and vp^o, respectively. Therefore 
the 4x4 Hamiltonian is given by 



0 ^ 

Hnik) ) 


( 1 ) 


with Hs{k) := uf [Xs(fc — A®) • cr — Ag], where s = 
L or i? is the index of the chirality and Xl{r) = 
and cr* is the Pauli matrix. Well, we calculate the cur¬ 
rent density in the presence of the external fields. The 
calculation is performed by the variation of the effective 
action, S'eff[Hi,e“], with respect to the gauge field, as 
j“(r) •=—{e'l{r)/\e{r)\){5Ses/5Ai{r)). The effective ac¬ 
tion is defined as 


g-Se„fo..e“] _ J p^p^texp(-5[V;,V'^A„e“]) , 

S[tlj,tlj\Ai,ei] ^ y drd^r {T,r)Cip{T,r) + c.c. 
£:=|e(r)|[a/5r-iJ(-iV„)], (2) 


where ip is the fermionic field, t and r denote the imag¬ 
inary time and spatial coordinate, respectively, and Ai 
is the vector potential. Here c.c. represents the complex 
conjugate combined with the change of the sign of the 
derivative operator djdr. Also, in Eq.(2), the Jacobian 
is given by |e(r)| := dete“(r), and the covariant deriva¬ 
tive is —iVa := ep{r){—idj — eAj(r)) with a = i. Using 
Eq.(2), we obtain that the current density up to the first 
order of the magnetic field and the TMF is given by 

reV(A^ - A^) „ , evpiX^ - A^)A ' 

^ ^ J ’ 

( 3 ) 

at zero temperature and up to the linear order in 
where the details of the calculations are described in 
Ref. [45]. Here, The vector representation of the TMF, 
T“, is defined by {T°‘)i := (l/2)e*-’'‘T“^. For the deriva¬ 
tion of Fq.(3), we introduced a momentum cutoff scheme 
I A; — A® I < A for the Weyl node of the chirality s. Phys¬ 
ically, A corresponds to the momentum range from the 
Weyl points in which the cone structures of the band of 
the lattice system is approved. 

The first term represents the CMF in the presence of 
the chiral chemical potential (i.e. Aq ^ Xq), and then 
reproduces the previous result for the CMF [46]. On the 
other hand, the second term in Fq.(3) is a new one, which 
raises the TCMF; i.e. the current is generated by the 
TMF for the pair of Weyl points which are shifted in the 



FIG. 1. (Color Online) Ground state current j induced by 
(a) edge and (b) screw dislocation with the Burgers vector b. 


momentum space due to broken time-reversal symmetry 
(TRS). This point is in sharp contrast to the usual CME, 
which requires breaking inversion symmetry. 

We comment on the relation between the TCME and 
the chiral anomaly. One may expect that when, Xjp = 
—A^, the TCME is described by the topological 0-term, 
which is the consequence of the chiral anomaly like the 
CME and anomalous Hall effect. However, there is no 
mixed chiral anomaly term of U{1) field strength and 
torsion in four-dimensional spacetime [40, 47]. This point 
is resolved by the observation that the Lorentz symmetry, 
which is postulated in the calculation scheme in Refs. [40, 
47], is broken in the cutoff scheme used for the derivation 
of the second term of Eq.(3), which is correctly applicable 
to realistic condensed matter systems. 

Now we discuss the consequences and physical pictures 
of the TCME. The TCME is realized in two types of 
lattice dislocations, (a) case of edge dislocation: oc 

AA^ Tf, and (b) case of screw dislocation: oc AXzTf, 

with AAo := Xpp—Xpp. Their schematic pictures are shown 
in FIG.l. These responses can be understood with the 
following semiclassical picture: Case (a): Edge disloca¬ 
tion is regarded as the (0,1,0)-“surface” of the extra lat¬ 
tice plane made up of the blue and green atoms in FIG.l- 
(a), which harbors a chiral Fermi arc, when two Weyl 
points are shifted in the fc^-direction. The electrons in 
the Fermi arc state are the very origin of the current in¬ 
duced by the edge dislocation. Case (b): There is a chiral 
Fermi arc mode on the dislocation line. The electrons in 
the mode rotate around the screw dislocation line, and 
due to the screw dislocation the rotating motion causes 
the current along the Burgers vector. 

The situation is similar to that of the three- 
dimensional integer quantum Hall state (3DIQHS) [48] 
with dislocation which is the staking of quantum Hall 
state layers characterized by the vector Gc = {2Trncf a)n, 
where Uc, a, and h are the first Chern number, the lat¬ 
tice constant, and the unit vector along the staking direc¬ 
tion. In the 3DIQHS, there are one dimensional n chiral 
modes along the dislocation line, when the topological 
number, n = bg ■ is nonzero [49]. This condition 

for the chiral modes is similar to that for the TCME, 

bg ■ (A^ - A«) ^ 0. 
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However, there are the following significant differences. 
The chiral modes of the 3DIQHS are exponentially lo¬ 
calized at the dislocation, and separated from the bulk 
higher-energy states, while, as will be shown below, the 
chiral modes of the WSMs exhibit power-law decay. We 
call them quasi-localized modes. Moreover their spec¬ 
trum is not isolated from the bulk spectrum but appears 
as its envelope as shown in FIG.2a, and therefore they 
can be easily mixed with the bulk modes. 

Spectral asymmetry and ground state current — We 
confirm the TCME due to dislocation by using an al¬ 
ternative approach other than the linear response the¬ 
ory based on (3). Our approach here is to calculate 
explicitly the spectrum and the eigenstates of the Weyl 
Hamiltonian with dislocation and the ground state cur¬ 
rent. We also show that the quasi-localized modes along 
the envelope of the bulk spectrum contribute to the ef¬ 
fect. For simplicity, we set vp = e = 1 and assume 
Xq = Xq = 0, and the Weyl points lie symmetrically on 
the fcz-axis, = Xz. In the presence of the 

screw dislocation at x = y = 0 along z-axis, of which 
Burgers vector is bg = —bgZ, the vielbeins are given by 
e| = —bgy/2TTp^, e? = bgxj^-np^, and e^ = for others, 
with p = + ?/2 [ 42 ]. Even with the dislocation, re¬ 

mains a good quantum number. Then, when is fixed, 
the Hamiltonian is equivalent to that of two dimensional 
massive Dirac model in the presence of the magnetic flux 
at the origin, whose amplitude is = kzbg, 


o-screw _ f rr± , _2:1 





(4) 


with the mass = kz—XsX. The equivalence of a screw 
dislocation and momentum-dependent magnetic field has 
also been pointed out in Refs. [40, 50]. 

The spectrum of consists of two types of eigen¬ 

states: one with the eigenenergies satisfying \E\ > \mf.J 
and the other one with E = irnl^^. The former does 
not contribute to the ground state current owing to the 
one-to-one correspondence between E+ > |m|^| and 
E~ < —\ml I modes as E+ = —E~, and between the 
states of Weyl nodes with the opposite chiralities [45]. 
On the other hand, the latter does contribute owing to 
asymmetry, i.e. the absence of one-to-one correspon¬ 
dence between E = ml^ and E = —ml^ modes. This 
asymmetry is called the parity anomaly [51, 52]. The 
asymmetric spectrum consists of discrete modes whose 
wavefunctions exhibit power-law decay, and continuum 
scattering modes which spread over the whole system 
[45, 53, 54]. The schematic picture of the density-of- 
state of the full spectrum is shown in FIG.2a. Moreover, 
the ground state current calculated from the asymmetric 
spectrum is = —Lz&gAA/27r^, which coincides with 
the expression obtained directly from the linear response 
theory (3) [45]. 



FIG. 2. (Color Online) (a) Schematic picture of the spectrum 
of the WSM with the screw dislocation. The black (white) 
bands along E = represent the relatively higher (lower) 

density-of-state compared with that of the opposite energy 
E = ■ (b) Lattice with a pair of screw dislocations with 

opposite Burgers vectors, (c) Numerical result for the spec¬ 
trum of the WSMs with a pair of screw dislocations. The blue 
curves are the envelope of the bulk spectrum. The opacity of 
the dots represents the expectation value of \p — Z'**® 

(see the gray scale bar). In the figures (cl-3), the modes with 
this value smaller than 0.1, 0.15, and 0.2 are plotted, (d) 
Current density along z-direction, jz{x,y). 


Numerical calculation — We confirm the spectrum 
asymmetry and the TCME for realistic lattice models by 
numerical calculations. We use the tight-binding model 
of WSMs [35] generalized to the case with dislocation. 


r 


E 


c 


t 

r+i+Si^zbg&(r) 


FT 


t—x,y,z 





r-\-i-\-6i^zbgQ{r) 


r*Cr. 




+ h.c., (5) 


where the 4 x 4-matrices, F*, satisfy the SO{5) Clif¬ 
ford algebra {FbF^} = 2S,j [55], F^ := \r\Vi]/2i, 
r = {x, y, z) and i denote the position of the atoms and 
the a;®-direction unit vector, respectively, and t,r, and 
d are the real parameters, and we suppose the lattice 
constant as 1 and lattice size x Ly x Lz- We intro¬ 
duced a pair of screw dislocations along z-direction with 
opposite Burgers vector at = ±(Z*®,0) as shown 

in FIG.2-b, by sliding the hopping directions in the first 
and third terms of Eq. (5) as 0(r) = — 1 for the re- 
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gion a; = 0, < y < while 0(r) = 0 for other 

regions. We numerically diagonalized this model and ob¬ 
tained the spectrums and current. Here the material pa¬ 
rameters are set as t = r = 1 and d = 3.6. The lattice 
constant is 1 and the amplitudes of the Burgers vectors is 
set as = 1. For the calculation, we imposed the open 
boundary condition along the x— and y—directions and 
periodic boundary condition along the z—direction, and 
set Lx = Ly = dZ*'* = 38 and = 100. 

As shown in FIG.2cl-3, we obtained the asymmet¬ 
ric spectrum in agreement with the analytic calculation. 
The asymmetric modes are localized at the dislocation 
line. The quasi-localized chiral modes are not isolated 
from the bulk but easily mixed with the bulk modes 
(FIG.2cl-3). The current density at zero temperature 
is shown in FIG.2d. We obtained the upward current 
along the screw dislocation and downward current along 
the anti-screw dislocation due to the TGME. The to¬ 
tal current per the unit length toward z-direction due 
to one dislocation line is J^/Tz = 0.087, which is cal¬ 
culated by the summation of the current density in the 
a: > 0 half-plane, and this value is in the same order 
as that estimated from the linear response theory (3), 
J^/Lz 0.1. For the estimation, we set the cutoff as 
A ^ l/(lattice constant) = 1. 

No-go theorem of CME — The existence of the TGME 
in the realistic lattice system may seem to contradict with 
the no-go theorems of the ground state current [35, 36]. 
However, they prohibit the total current, but not the 
local current density. Therefore, the current along the 
dislocation line can exist, as we found[45]. 

Experimental implication — Here we present two ex¬ 
perimental setups to observe the TGME in TRS-broken 
WSMs, for which Eu 2 lr 207 [22] and YbMnBi2[23] are 
candidate materials. The first one is a scanning SQUID 
measurement, which can detect weak inhomogeneous 
magnetic fields [56, 57]. If there is a pair of disloca¬ 
tion, the circulating current occurs. The magnitudes of 
the current and the induced magnetic field are estimated 
as / ^ 10“®A and B ^ 10”^T, respectively, for both 
Eu 2 lr 207 and YbMnBi 2 . Here we used Eq.(3) and the 
material parameters, vp ^ lO^m/s and a ^ lOA and set 
bg = a and A ^ A ^ 1/a, where a is the lattice con¬ 
stant. Also, for the estimation of the magnitude of the 
magnetic field, we used a typical value of inter-distance 
between dislocations, 10®A [58]. It is feasible to detect 
B ^ 10”^T via the scanning SQUID. 

The second one is a nonlocal transport phenomenon, 
which was observed in quantum Hall materials [59, 60]. 
The experimental setup is shown in FIG.3. If the bulk 
contributions are completely negligible and there are 
only the chiral modes at the dislocation lines, V 34 := 
V 3 — V 4 = 0 despite /12 A Oj then, the nonlocal resistiv¬ 
ity i?i 2,34 := V 34 /I 12 is equal to zero [59]. On the other 
hand, if the nonlocal transport is negligible, V 34 > V 3 ' 4 ' 
holds for L 1 L 3 < L 1 L 31 when I 12 A 0- Therefore, if 


V3, 
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FIG. 3. (Color Online) Experimental setup for the nonlo¬ 
cal transport due to the TGME. The thick red solid (blue 
dashed) line represents the (anti-)dislocation line. The leads 
are attached at the black points, Li{i = 1, 2, 3, 4, 3/4'). The 
line L 4 L 2 , L 3 L 4 , and L 3 /L 4 / are parallel and have the same 
length. Here we suppose L 1 L 3 < L 4 L 31 . Vi is the voltage at 
Li, and I 12 is the current. 

.Ri 2,34 < 7 ?i 2 , 3 ' 4 ' is observed, it is the fingerprint of 
the chiral current due to the TGME. The effect can be 
discriminated from any previously reported conventional 
transport induced by dislocation[58, 61-69]. 

We also comment on effects of impurities and disori¬ 
entation of the dislocation. First, the current due to the 
TGME is expected to be robust against weak disorder. 
It is because that Eq.(3) is independent of the scatter¬ 
ing time, like the intrinsic contribution to the anomalous 
Hall effect [70]. More precisely, the current is due to the 
edge modes in the Fermi arc on the surface of WSMs, and 
these modes are supported by the Weyl points, which are 
protected by the Chern number, and hence, robust againt 
weak disorder. 

Next, in real experimental setups, it is difficult to align 
the dislocation line orthogonal (parallel) to the line con¬ 
necting the Weyl nodes exactly in the case (a) (case (b)). 
Even when they are not orthogonal (parallel), as long as 
they are not parallel (orthogonal), the current parallel to 
the dislocation line still exists. Supposing that the dislo¬ 
cation line is parallel to the z-axis, the current is given 
by = evpL^AiX^ - A^) • 6 g/ 47 r 2 in the both cases (a) 
and (b). 

Summary — In this letter, we have discussed the 
TGME in WSMs caused by dislocation. We have con¬ 
firmed that it is possible to occur and experimentally 
observable in realistic materials, and argued that the 
Lorentz symmetry breaking is important for it. 
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Supplemental Material 


Derivation of Eq. (3) 

In this section, we derive the expression for the current 
in the presence of the magnetic and torsional magnetic 
responses of the current density, Eq.(3). The derivation 
consists of two steps: first we derive the expression for 
the Green function in the presence of the gauge field and 
vielbein using the gradient expansion Eq.(S-4), and next 
we calculate the current density by using Eq.(S-4) and 
obtain Eq.(S-18), which is equivalent to Eq.(3). 

First, we calculate the single-electron Green function. 
The Green function in the presence of the gauge field and 
vielbein which is defined by 

G'(ri,ri,T 2 ,r 2 ) 

_ / T)'ipT)'ip'ftp{Ti,ri)'il^^T 2 ,r 2 )ex.p e“]) 

exp{-S[ip,'tp'i,Ai,ef]) 

(S-1) 

Then, the following differential equation holds : 
i C{eN,ri,-idri)G{eN,ri,r2) 

+ G{eN,ri,r2)t*{-eN,r2,-idr2) = - r2), 


with C(eN,r,-idr) := |e(r)|[ieAr - i?(-iVa)] - eAo], 
and |e(r)| := dete“(r) Here Sn = (2fV -|- l)7rT is the 
Fermionic Matsubara frequency with the temperature 
T, and G{sN,ri,r2) := Jq G(t, ri,0,r 2 )e“''^'^'^dr is the 
Fourier component of the Green function. Now, using 
the spatial Wigner transformation defined as f{R,p) := 
f d^re~‘'^'Pf(R + rl2,R— r/2), Eq. (S-2) is rewritten 
into 

+G(eiv, H, p)] = 1. (S-3) 

In the gradient expansion up to the first order in diAj or 
9ie“, the Green function becomes 

G{eN,R,p) = &°\eN,R,p) + &^\eN,R.p) + --- , 


&°\eN,R,p) = ’ 

&^\eN,R.p) 

i (9£o(eAf,7r) 1 dCo{eN,T^) 

■ Gq {en, tv) [eFab{R) + T^i,{R)tTc] L^=e^(ii)(p._eAi(R)) ’ 

(S-4) 


Here TVa := e\{R){pi — eAi{R)) is the gauge-invariant 
mechanical momentum, while pi is the canonical momen¬ 
tum, and the field strength and the torsion with the in¬ 
dices of the local orthogonal coordinate are, respectively, 
defined as Fab '■= eaSlFij and := e\e{Tfj. The free 
Lagrangian is defined as Co{eN,Tv) := ien — 

Next, using Eq. (S-4), we calculate the current density 
and derive Eq.(2). The current density are defined by 
fir) :=-{eflr))/\e(r))\){5Ses/SAi{r)). Therefore 


j\R) 


J'D'ilj'D^p'‘^p^T,ri)^-^tp{T,r2)e ^ 


f 'Dip'D'ip'^ e~^ 




= —VTr 
2 ^ 

N 

G{eN,r2,ri) 

Me/; 


dCo{eN,Tv) 

divi 

,T’ 2 — ^ 


— —eAi{r 2 )) 


d^P ^ r 5Eo(ejv , tt) 
(27r)^ _ diri 


\7Ta^ei^{R){pi-eAi{R)) 


Ai^^-K^)G{sN,R,p) 


Here Tr means the trace over the band indices, we used 
that dtjbAi = e\e{r)\e\{r)dLQ/diVa and the second line 
does not depend on t due to imaginary time-translation 
symmetry. Note that tv of the third line is the operator 
though that of the fourth one is the c-number. Using 
Eq.(S-4), up to the first order in diAj or diCj, the ex¬ 
pression of the current density becomes 

(S-6) 

with the zeroth-order terms. 


eT 


dCo{eN,Tv) 

dwi 


&'^\eN,R,p) + c.c. 
and the first-order terms, 

= V / ^^Tr [ - 

2 (27r)3 

G«(£^,i^,p)l 


TTa — Ga (■^)) 


dCo{eN,Tv) 

dni 


'^a — e'^iR) {Pi—eAi{R)) 


dCo{eN,Tv) 

d fi'n-vr 


\7Ta—e'^^{R) {pi — eAi{R)) 


-%^)&°\£n,R,p) +C.C.. 

The zeroth-order terms (S-7) can be rewritten as 
f d^TV , 


P^°'>{R) = eJ 2 


^v'^’^{Tv)nFi£n,7v)- 
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For the derivation, we inserted the iden¬ 
tity, It, = between dCo/diri 

and G(o) in Eq.(S-7), and used the formula, 
E“=-oo [l/(ii^JV -t) + l/(-i£jv - t)] = (1 - 2nF{t))/T, 
for the summation over the Matsubara frequency, 
and f (Pp = |e(fi)| f (Pn. Here, n is the band index, 
is the energy, u”’“( 7 r) := den,TvldTTa is the group velocity, 
ni?(e) := 1/(6^/"^ -I- 1) is the Fermi distribution function, 
and |m^) is the Bloch state. This term corresponds to 
the summation of all contributions to the current from 
the electrons in the occupied states in the absence of 
magnetic and torsional magnetic field. 

Now, we move on the calculation of Eq.(S- 8 ). The 
sum of the second term of Eq.(S- 8 ) and its complex con¬ 
jugate is zero, since [•••]* = [• • • Then, by 

using Eq.(S-4), we obtain 


- 


E 

N ' 


£0 ^(£JV,7r) 


-I- C.C.. 


' (27r)3 ’’ 

9£o(£jv,7T 

diTh 


(9£o(£Ar,7r) 


dwi 

-£o ^(£Ar,7r) 


„-i, a£o(£Ar,7r 

£0 

(S-10) 


Moreover, inserting the identities, 1,, = |u”) (u^|, we 

obtain 


-ieT 


E 


X ( I 


N,n,m 

dH 

n 


d^TT 


diTb 


dH 

dwi 

1 


m 


m 


dH 

dlTa 


(i£Ar - £n)(ieAf - £m)(ieAr “ £l) 


X [eFab{R) + T^b{R)T^c] + C.C., 


(S-11) 


where the indices tt are omitted like £„ := £n,-K and 
\n) := \u'^). There are three types of contributions to 
the summation over the band indices n,m,h (a) all the 
three are the same, (b) two of them are the same and 
the other is different, and (c) each one is different respec¬ 
tively. However the contribution (a) is found to be zero 
because of the antisymmetry of [eFab{R) +T^b{^)'’^c\ 
under a ^ b. Moreover, the contribution (c) is also zero, 
since our model of the WSM, Eq.(l), consists of two two- 
band Hamiltonians independent of each other, and then 
the overlap of three or more bands is zero. Therefore, 
we have only to consider the contribution (b), and then 
obtain 

• p >-73 

= E J (27r)3 + ^abi + -^61a) 

(S-12) 


with 


Mabc ■■= v^npien) {n, & |(£„ - i7) | n, c) 

-I- u” (n, b I npiH) | n, c) - v”n_F(£„) (n, & | n, c), 

(S-13) 


where we used the abridged notation, |n,a) := 

For the derivation of Eqs.(S-12,S-13) we used the formu- 

] Y^oo _1_ _ tnp{s)—snp{t)-\-nF{s)—nF{t) 

2^N— — oo (isN—ty^{'iSN—s) T{t—sy^ ’ 

(n I dH/diTa I m) = (£’" — £”) (n | m, a) for n m, and 
J2mfi^m){n,b\n,c) = {n,b \ f{H)\n,c) for any func¬ 
tion /. Moreover, using the relationship (Mabc)* = Macb, 
we obtain 

= ^ E / ^ ["^23 + Tin,] S^^^M^rbc, (S-14) 

where is the antisymmetric symbol. Furthermore, 
since 'u"n^(£n) = dupisn)!dna and only antisymmetric 
parts of Mabc contribute, using integration by parts, we 
find 

= ie ^ y [eF 23 + T^nd] £“^°u>F(e„) (n, 61 n, c) 

+ fE/ -^^T^'iS''^’^SnnF{£n){nMn,c). (S-15) 

Using the Berry curvature is defined by := 

—i£“^'^ (n, &|n, c), and the vector representation of the 
TMF, r“, is defined by T“ := it can be 

rewritten as 

jTF(^r) 

= -eE / (e-Bi+£i“7ra)ni.(£„) 

- i E / ^Kner^UFiSa). (S-16) 

It is noted that the term containing Bi is equal to 
the expression for the CME derived by Son and Ya¬ 
mamoto [46], and the others are new terms that rep¬ 
resent the current induced by the torsion. Neglecting 
the last term, which is, as we will discuss later, less im¬ 
portant than the others in the case of WSMs, Eq.(S- 
16) can also shortly derived from the substitution of 
the magnetic field or the field strength in the absence 
of the vielbein, —i[(—ii92 — e>42),(—i^a — eA^)] = eBi 
with the field strength in the presence of the vielbein, 
-i[-iV 2 ,-iV 3 ] = eBj -f Tj“(-iVa), where [U,V] := 
UV — VU is the commutator. This justifies the anal¬ 
ogy between the TMF and the magnetic held. 


-I- C.C., 
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Finally, we substitute the energy, group velocity, and 
Berry curvature of the model of the WSM, (1), into 
Eq.(S-16) and derive Eq.(3). We characterize the four 
bands of the Hamiltonian (1) as n = (s,±), with s = 
L or R, where s is the index of the chirality and +(—) 
means the higher (lower) band of the Weyl cone. Then, 
their energy, group velocity, and Berry curvature are 
given by 


e‘^’^{k) = VF [±|fc-A®| - Ag] 


,S,± 


(fe)=± 


VFik-X^) 


= ±x 


fc-A«| 
k-X^ 
’ 2 |fe-A^|3’ 


(S-17) 


Using Eqs.(S-16, S-17), we obtain 


Ag') 

47r^ 


Bi{R) 


evFjX^- X ^)A 
Att'^ 


Tf(H) 


(S-18) 


at zero temperature and up to the linear order in 
For the derivation of Eq.(S-18), we introduced a momen¬ 
tum cutoff scheme \k — X^\ < A for the Weyl node of 
the chirality s. Physically, A corresponds to the momen¬ 
tum range from the Weyl points in which the cone struc¬ 
tures of the band of the lattice system is approved. Note 
that the last term of Eq.(S-16) yields second(or more)- 
order contributions in , and then less important as 
mentioned before. Eq.(S-18) is the correction of current 
due to the TMF and magnetic field and is equivalent to 
Eq.(3), then the derivation of Eq.(3) has been completed. 


Ground state current in the presence of screw 
dislocation: analytical calculation 


In this section, we calculate the ground state current 
raised by the TCME in the case of screw dislocation, by 
calculating directly the eigenstates of the Hamiltonian 
with the torsion. This is an alternative approach for the 
derivation of the TCME, which does not rely on Eq. (3). 
For this purpose, we, first, analyze the spectrum of the 
Hamiltonian (Eq.(4) in the main text). 


rrscrew 

Ht = 


= Xs [Ht +^, 0 -^] : 





(S-19) 


where - XsA, = k^bg, p = y/x'^ + 

Xl{r) = and cr* is the Pauli matrix. 

For the calculation of the spectrum, it is useful to clar¬ 
ify the symmetry of the eigenstates of 77^. Suppose 


|k)j, the eigenstate of with eigenvalue k. Since 
{i7£,cr^} = 0, where {U,V} := UV + VU is the an¬ 
ticommutator,the state |k)^ is also the eigenstate 
with eigenvalue —n. Therefore, we can choose the eigen¬ 
functions to preserve the doublet structure, |—k)^. = 

cr* |k)^ , for k 0. On the other hand, there is no double 
structure in the zero eigenstates. Since the hermitian op¬ 
erator maps zero eigenstates of to zero eigenstates 
of , then we can choose the zero eigenstates also as 
eigenstates of , denoted by with = 

<^i and Ui = ±1. There is another symmetrical 

property between the eigenstates of 77^ with different 
kz- Since the transformation kz —>■ —kz corresponds to 
the flip of the direction of the effective magnetic field, 
077^0“^ = holds, where 0 = ia^K is the time- 

reversal operator for spin-1/2 fermions and K is the com¬ 
plex conjugation operator [71]. Therefore we can impose 
\k)_^ = Q\k)^^ and = 010*,^^)^ , because 

of{cr^0} = O. " 


The eigenstates of 77|//®'^ can be constructed 
from and jOcrJ^,^. Indeed, := 

+ c£i^(K) |-K)fe^, with K > 0, and 

IOctJj. are the full spectrum of 77|“®/^, with eigen¬ 
values and respectively. Here 

the coefficients are given by c^’^(k)) = 

(4(«2 + (mt)^))-^/^(±sgn(m^J((«2 + (m|j2)V2 ± 
k)^/^, ((K^-|-(m(j^)^)^/^=FK)^/^)). Moreover, |V'£’^(«^)) := 
0 |V'-/^(^)) and are the eigenstates of 

with eigenvalues )‘^ and respec¬ 

tively. 


Now, we calculate the ground state current in the pres¬ 
ence of the screw dislocation. As yet, we have not distin¬ 
guished discrete and continuum states. From now on, we 
use Ki to express the discrete eigenvalues of 77/- and (k, /) 
to label the continuum states, where k is the continuum 
energy eigenvalue, and I is a discrete quantum number, 
e.g., the angular momentum. The current operator is 
defined by 977-/7/97:, = x«a^ + Xs{-ibgy/2Trp^)a- + 
{bgx/2TTp‘^)ay}. At least up to the first order in bg, the 
correction to the current operator due to the dislocation, 
i.e. the second and third term above, does not contribute 
to the expectation value because these terms are odd un¬ 
der the transformation x —>■ —x or y —>■ —y. Therefore, 
the current is the sum of the expectation values of XsO’^ 
for the occupied states which consist of discrete nonzero, 



discrete zero, and continuum states, and then we obtain. 


Absence of total current and possibility of local 
current 




L,dkz 




+ X! (0*.<TjXsCr^|0*.ai) 




i-.a-iXs-mf,^ <0 
^OO 


pOO _ 

Jo , 


L,dk, 






/ oo 

dK^{K,l\a^\K,l)^ 

-OO I 


(S-20) 


where Lz is the size of the system. Here we introduce the 
momentum cutoff scheme, |m^ | < A, i.e. the domain of 
the integration is the same as that used in the calculation 
of Eq.(4). The first term in the square braket is equal 
to zero, since orthogonal to 

The second term is the index of the Dirac operator, , 
which is an integer and the difference in the number of 
its normalizable zero modes with cr^ = +1 and = 
— 1. The index is given by Nk^ := —sgn($fc^)[|$fc^|/27rj 
[53, 54]. The normalizable zero modes exhibit power-law 
decay for large distance from the dislocation; i.e. they 
behave like |0q_i) oc —iy)*"^) for > 0, 

and |0i^+i) oc {x + iyy~^,Q) for 4)^^ < 0, where 

i = 1,2,-•• ,|iVfc^| [53, 54]. Now, we move on to the 
third term of Eq.(S-20). One may expect that it is equal 
to zero, since [k, 1 ) = ]—k, 1 ), is orthogonal to [k, 1 ) for 
almost all values of k. However, the scattering states near 
K = 0 (their amplitudes oc cos{kp + Si)/-^ with 6 i the 
phase shift) cause a delta function peak of {K,l\a^\K, l)j.^ 
at K = 0. Indeed, from an explicit calculation [53], it 
has been shown that {K,l\a^\K,l )= Ckzd{n), with 
Ck^ = 4’fc^/27r — Nk ^, and then the third term is equal to 
Cfc^. Substituting them into Eq.(S-20), we obtain 




Lzdkz <^k 

-A-A 27 r 


27r 


LzbqAX 

27r2 


(S-21) 


which is coincident with the expression obtained directly 
from Eq.(3) by the following reason. In the presence 
of the screw dislocation with the Burgers vector —bgZ 
the torsion is given by T| = T^y = bg5^'^\x,y). There¬ 
fore, the total current derived from Eq.(3) is = 
—LzevpiXf — X^)Abg/4Tr‘^. In this section we have set 
X^ = —Xf = X and e = vp = 1, and therefore we obtain 
= —LzXAbg/2'K'^, which reproduces Eq.(S-21). 


In this section, we show that, if the system is periodic 
in a certain direction, the total current is always zero, 
even in the presence of a magnetic field or lattice stain 
and dislocations, while the local current is not. The argu¬ 
ment is the extension of that presented in Ref. [35]. We 
start with the general Hamiltonian of electrons in solids 
(set e = I in this section): 

H = I d^r^ (-iVi - Mp, z) f + V{p, z), (S-22) 

where V is the potential term, in which the effect of the 
dislocation is included. Here i = x,y,z, r = {p,z), 
and p = (x,y), and we impose the periodicity in the 
z-direction: 


Ai{p,z) = Ai{p,z + a),V{p,z) = V{p,z + a). (S-23) 


Suppose tp is one eigenfunction of the Hamiltonian 
and define the Bloch wave function ipn,k^ {p, z) = 
(p, z), whose energy is en,fe^. The total current 
along the z-direction is given by 

= j '^n,kSP,z)^yn,kAP,z)nF{Sn 


dk 


dHu 


= -Yl I Ki,kAP,z)-^Un,kAP,AnF{en. 


! BZ 


dkz 


(S-24) 


where Hk^ = e and uf is the Eermi distribu¬ 

tion function. Here we use the identity : 

J d^r <,fe^(p,2:)^^u„,fe,(p,z) 

■ J d^r ul k, (p, z)Hk,Un,k, (p, z) 


A. 

dkz 
dSn^kz 


dkz ’ 

which follows from 


d 

dkz 


J d^r uAkziP^ z)un,kziP,z) 


A. 

dk. 


(S-25) 


I = 0, (S-26) 


and then we can rewrite Eq. (S-24) into 
dkz dCji kz 


j. = -E 


f BZ 


27r dkz 


~Tlp(^£n.kz ) 


=-si: i: 




ds npis) (S-27) 




^i -1 


Here for each region k € {kAi^^A)^ monotonically 
increases or decreases, and AA = 0 and = 27r/a. We 
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find that Eq.(S-27) is always equal to zero owing to the 
periodicity of the dispersion in the wave number space, 
^n,k=o = ^n,k= 2 Tr/a- Therefore, we found that the to¬ 
tal current along the z-direction is zero. In the above 
derivation, it is essential that the integrand with respect 
to kz can be rewritten into the total derivative with re¬ 
spect to kz, and this key factor follows from the fact that 
the integral over the real space of \un,kApT ^)\'^ is equal 
to 1, which results in Eq. (S-26). Instead, without the 
integration over the real space, 

^ K,fc, (P, (P, z)] ^ 0. (S-28) 

Then, in the case of local current, the above argument in 
the case of the total current does not hold. Hence, the 
local current is not always zero, unlike the total current. 
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